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measure for the metal mass? 

W. Kapferer^, T. Kronberger^, J. Weratschnig^ , and S. Schindler^ 



o 
o 

(N 



6 



> 

cn 

o 
i> 
o 



X 



Institut fiir Astro- und Teilchenphysik, Universitat Innsbruck, Technikerstrasse 25, A-6020 Innsbruck 

ABSTRACT 

Aims. We investigate whether X-ray observations map heavy elements in the Intra-Cluster Medium (ICM) well and 
whether the X-ray observations yield good estimates for the metal mass, with respect to predictions on transport mech- 
anisms of heavy elements from galaxies into the ICM. We further test the accuracy of simulated metallicity maps. 
Methods. We extract synthetic X-ray spectra from N-body/hydrodynamic simulations including metal enrichment pro- 
cesses, which we then analyse with the same methods as are applied to observations. By changing the metal distribution 
in the simulated galaxy clusters, we investigate the dependence of the overall metallicity as a function of the metal 
distribution. In addition we investigate the difference of X-ray weighted metal maps produced by simulations and metal 
maps extracted from artificial X-ray spectra, which we calculate with SPEX2.0 and analyse with XSPEC12.0. 
Results. The overall metallicity depends strongly on the distribution of metals within the galaxy cluster. The more 
inhomogeneously the metals are distributed within the cluster, the less accurate is the metallicity as a measure for 
the true metal mass. The true metal mass is generally underestimated by X-ray observations. The difference between 
the X-ray weighted metal maps and the metal maps from synthetic X-ray spectra is on average less than 7% in the 
temperature regime above T > 3 x 10^ K, i.e. X-ray weighted metal maps can be well used for comparison with observed 
metal maps. Extracting the metal mass in the central parts (r<500 kpc) of galaxy clusters with X-ray observations 
results in metal mass underestimates up to a factor of three. 
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1. Introduction 

Since the first observations of the 7 keV iron hne feature 
in the 1970's by Mitchell et al. (1976) it has been evident 
that the intra-cluster medium (ICM) does contain gas al- 
ready processed by stars. X-ray spectra are the only mea- 
sure for the metallicity of the ICM. With X-ray observato- 
ries like XMM-Newton or Chandra it is nowadays possible 
to extract metallicities in certain regions of a galaxy clus- 
ter and construct metallicity profiles and X-ray weighted 
metallicity maps (Pratt et al. 2006, Durrett et al. 2005). In 
addition to observations galaxy cluster simulations includ- 
ing gas and dark matter (DM) physics are an ideal tool 
to investigate the evolution of these systems. By includ- 
ing different enrichment processes it is possible to trace the 
metals in the ICM. Recent work in this field (Cora 2006, 
Kapferer et al. 2006, Kapferer et al. 2007, Domainko et 
al. 2006, Moll et al. 2007, Springel & Hernquist 2003) are 
able to reproduce the inhomogeneous distribution of met- 
als within the ICM as observed (Durret et al. 2005). Given 
the large mass fraction of the ICM in a cluster (15-20%) 
compared to the mass fraction of the galaxies (3-5%) a lot 
of metals must have been produced within the galaxies and 
then been transported together with part of the interstellar 
medium (ISM) into the ICM. X-ray observations can also 
distinguish the lines of different elements, e.g. elements like 
Si, S, O from core collapse supernovae or the elements Fe 
and Ni from supernovae la (Baumgartner et al. 2005, Ettori 
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et al. 2002, Sanders et al. 2004, Finoguenov et al. 2002) and 
can hence give information on the origin of the metals and 
interaction in clusters. From these observations we know 
(e.g. Mushotzky 1999), that most of the metals are in the 
ICM and not in stars. Thus, to trace metal formation and 
evolution it is crucial to get correct estimates for the metal 
content in the ICM. 

Three dimensional hydrodynamic galaxy cluster simu- 
lations give us an estimate of the 'true' metallicity distribu- 
tion within simulated galaxy clusters. By constructing arti- 
ficial X-ray profiles and maps it is possible to test, whether 
metallicities obtained by X-ray spectra are a good measure 
for the real metal content in the ICM or not. 

The paper is structured as follows. First we compare 
X-ray weighted metal maps with metal maps obtained 
from synthetic spectra, which we construct with SPEX 2.0 
(Kaastra et al. 1996) and investigate with XSPEC 12.0 
(Dorman et al. 2003). Then we compare different heavy 
element distributions in the ICM and their spectroscopic 
appearance. We give constraints on the accuracy of spec- 
troscopic metallicities as a measure of the real metal content 
in the ICM. Throughout the paper we assume solar ratios 
(Anders & Grevesse, 1989). 



2. The Model Clusters 

We take two realistic ICM distributions from model clusters 
for all further calculations into account. We apply a stan- 
dard ACDM simulation with gas physics and constraints 
random fields (see Kapferer et al. 2007, Kapferer et al. 2006, 
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Domainko et al. 2006 and Schindler et al. 2005 for more de- 
tails) to obtain a reasonable ICM gas distribution. The first 
model cluster is a merging galaxy cluster, hereafter Model 
Cluster A, whereas the second system is a more massive 
galaxy cluster, hereafter Model Cluster B. The properties 
of the model galaxy clusters are as follows: 

— Model Cluster A: The cluster forms at z^l.5 and has 
two major merger events at z=0.8 and z=0.5. The final 
total mass is 1.5 x lO^** M© within a sphere of radius 1 
Mpc. 

— Model Cluster B: The formation redshift for this clus- 
ter is z^l.7. It shows four minor merger events a z=1.4, 
z=l.l, z=0.5 and z=0.3. The cluster has a final mass of 
3.4x10^'* Mq within a sphere of radius 1 Mpc. 

The galaxy cluster simulations have two enrichment 
processes included, namely galactic winds and ram-pressure 
stripping. The enrichment process description reproduces 
the metallicity profile of non-cooling flow and massive 
galaxy clusters (see Kapferer et al. (2007) for more details 
about the enrichment history and a detailed analysis on 
the metal distribution) . With these two realistic models for 
ICM distribution we performed all analysis on the differ- 
ence of the true metal mass, which we know exactly from 
the simulation, and the metal mass we would obtain from 
X-ray observations of our model clusters. 

3. How to construct X-ray weighted metallicity 
maps 

So far X-ray weighted metal maps where used to obtain 
metal maps from simulations. We want to test here how 
accurate these maps arc compared to maps, that are ob- 
tained by using spectra. The total emissivity by thermal 
Bremsstrahlung for the ICM above T > 3 x 10^ K is ap- 
proximately given by 
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where T is the ICM gas temperature, rig, n.i are the number 
densities of electrons and ions, Z is the number of protons of 
the bending charge and gs the frequency averaged Gaunt 
factor. The emissivity is mainly dominated by the density of 
the ICM. Besides the continuum emission several emission 
lines in the ICM can be observed. The most prominent lines 
in the ICM are typically in the 7 keV iron line complex. 
These lines origin from many stages of ionisation of iron. 
The power emitted by the most prominent emission lines 
in the ICM, the iron K line complex, which is a resonance 
line, is given by 



Pjk = 8.63 X 10*^ {nenhA,iX,{t)Ejk) x 
Stt 
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where A^i is the abundance of the element el relative to 
hydrogen and Xi the fraction of the ion i, Ejk the transition 
energy from level j to fc, Egj the transition energy from 
ground level g to \, the averaged Gaunt factor g{T) and 
f'gj the effective oscillator strength. For more details see 
Sutherland & Dopita (1993) and references therein. Again 



the dependence of the emission on the density of the ICM 
and the temperature of the gas is evident. For a given set 
of fixed ICM density and iron abundance the temperature 
dependence of the emission ratio of a given iron line and 
the bremsstrahlung continuum can be written as 

«m.«i.n22i4ZI) (3) 

where VL is an arbitrary normalisation constant and A = 
Egj/k. Please note that this approximation is different for 
each single line of a given ion and that the temperature 
dependence of ^^(r), g{T) and 5s(T) is ignored here. 
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Fig. 1. The temperature dependency of Eq. 3, i.e. the 
ratio of the emissivity for a given line and the thermal 
bremsstrahlung continuum in the temperature range 0.3 
- 2x10^ K. 



In Fig. [T] the dependence of a given line and thermal 
bremsstrahlung emission as a function of temperature is 
given. In our model setup we enrich a fixed space with 
a fixed mass of metal. Then we construct X-ray emission 
weighted metallicities, see Eq 3. This is done by integrating 
the metallicities along the line of sight and the weighting 
by the square of the density and the derived temperature 
dependence, i.e. 



Metallicityx^r 



^ EiosnLMiosR(T)io, 

ay weighted 2 r)/'T~\, 



(4) 



where njos, M/o^ and R(T)ios are the density, metallicity 
and temperature-dependent weighting factor in each cell 
along the line of sight. By summing also perpendicular to 
the line of sight, we can obtain the mean metallicity over 
the whole cluster. 



4. Spectroscopic metal maps and the comparison 
with X-ray weighted metal maps. 

In order to simulate the X-ray emissivity of our model clus- 
ter we use SPEX 2.0 to construct thermal bremsstrahlung 
spectra including line emission. In Fig. [3] a spectrum for the 
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Fig. 2. Left: X-ray weighted metal map of the model galaxy 
cluster at z=0, the metallicity in each cell is weighted by 
the X-ray emission. Right: Synthetic metal map, using a 
spectrum for each cell. The image has 1.5 Mpc h^^ on the 
side. 
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Fig. 4. Ratio of X-ray weighted metal map to synthetic 
spectroscopic metal map. The image has 1.5 Mpc h"-'^ on 
the side. Histogram of the left image. The X-ray weighted 
metal map differs on average by 7% from the spectroscopic 
metal map. 



whole model cluster is given. For each cell of our simula- 
tion, which has a resolution of 128'^, we calculate a model 
spectrum which we than integrate along the line of sight. 
By using XSPEC we fit a spectrum to each pixel and hence 
construct a metallicity map, as shown in Fig. [2j 
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Fig. 3. Model spectrum of our model cluster, obtained with 
SPEX2.0. The count rate of the whole cluster is shown. 



The model we used in XSPEC12.0 is exactly the same 
which we used for constructing the synthetic spectra. In 
order to investigate the quality of our X-ray weighted 
metal maps we divided the two metal maps, the metal 
map obtained with XSPEC12.0 by fitting for each pixel 
the synthetic spectra with a MEKAL-model (thermal 
bremsstrahlung with line radiation for an optically thin 
plasma) (Mewe, Gronenschild & van den Oord, 1985) and 
the X-ray weighted metal map. The result is shown in 
Fig. 131 beside the division image a histogram of the ra- 
tios is given. The difference of the two maps is on average 
7%, leading to the conclusion, that X-ray weighted metal 
maps are a good approximation for spectroscopic metallic- 
ity maps, despite all the approximation made. In the out- 
skirts where the density and temperature of the plasma is 
lower, the deviation is larger, due to lower metallicities in 
the X-ray weighted metal. This can be understood in terms 
of the validity of the thermal bremsstrahlung approxima- 
tion (Eq. 1), which is only valid for a plasma with tem- 



peratures above 3 x 10^ K. In the outskirts of our model 
cluster r>600 kpc the temperature can drop below 3 x 10^ 
K, resulting in a too low metallicity. In this case the syn- 
thetic spectra give better results, because the Mekal model 
is valid for plasma temperatures above 3 x lO'' K (Mewe, 
Gronenschild & van den Oord, 1985). In regions 1 Mpc 
around massive cluster centres the temperature does typi- 
cally not drop below 3 x 10'' K, therefore the X-ray weighted 
metal maps used so far to obtain metal maps from simu- 
lated clusters are a good approximation. 

5. Results 

5.1. The mean metallicity as a function of the metal 
distribution 

As metallicities in the ICM are always obtained by aver- 
aging over areas or annuli in the field of view, the ques- 
tion arises, how the real distribution of metals infiuences 
the mean metallicity of an observed cluster. Simulations 
are the only way to address this question, because the 3D 
metal distribution is exactly known. In order to investi- 
gate how sensitive the mean metallicity of the ICM acts 
on the distribution of heavy elements we investigate an ex- 
tremely inhomogeneous metal distribution in the ICM. We 
add metals in a region with ~ 120 kpc on a side containing 
1.17x10* Mq heavy elements, which we place at different 
positions along a line perpendicular to the line of sight, 
which passes through the cluster centre, i.e. we always add 
the same amount of metals to the ICM and see how the 
mean metallicity of the galaxy cluster in an 1 Mpc radius 
area depends on the position of the metals. Note that the 
ICM density distribution corresponds to model cluster A, 
with zero metallicity, except for the metal box mentioned 
above. In Fig. [H] the projected mean metallicity of the ICM 
in an area with 1 Mpc radius around the cluster centre is 
shown as a function of the position of the cube with respect 
to the cluster centre. The ratio between the metal mass ob- 
tained from the mean metallicity and the true metal mass 
is presented as well. The geometry of the observation and 
the metal cube position is shown in Fig. [S] The metallicity 
measurements always refer to observations along the line of 
sight. 

It is evident that (see Fig. [H]) the position of the cube 
infiuences the result dramatically. This can be understood 
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Fig. 5. The geometry of the observation and the metal blob 
positions. The metal blob is placed at different positions 
along the plotted line, which is perpendicular to the line of 
sight. 
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Fig. 6. Overall metallicity of the ICM in a circle with 1 Mpc 
radius around the cluster centre as a function of the position 
of a region containing heavy elements. A metal blob with 
117 kpc on a side containing 1.17x10* Mq heavy elements 
moves through the cluster centre, as highlighted in Fig. [5l 
On the right side the ratio of the metal mass obtained from 
the mean metallicity and the true 3D metal mass is given. 
The true metal mass is always higher then the metal mass 
obtained from X-ray observations. Note that the cluster has 
zero metallicity, expect for the metal box . 



by the density distribution within the ICM. The line flux 
of a given element can be written as 

Fune J f{T)n,ndldA ^ J f{T)n,dM,, (5) 

where rie, rii are the number densities of electrons and 
the given element ions, /(T) is a temperature dependent 
emissivity and dMi is the differential mass of the given el- 
ement. In the centre the high density of the ICM results 
in the highest mean metallicity, which contributes most to 
the X-ray emission. When the metal blob is placed at a ra- 
dius of 1 Mpc the mean metallicity drops nearly one order 
of magnitude, because the outskirts do contribute less to 
the whole cluster emission. From this case of a metal blob 
placed in primordial ICM we can conclude, that the mean 



distribution depends strongly on the position of the metal 
blob within the ICM. 

Recalculating the mass of metals residing in the ICM from 
the mean metallicity, shows the same dependence, see Fig. 
[51 The mean metallicity always underestimated the metal 
mass present in the model cluster. The discrepancy is a fac- 
tor of 10 in the centre and increases to a factor of 100 in the 
outskirts. The reason therefore is again the fact that emis- 
sivity by thermal bremsstrahlung is mainly driven by the 
density of the ICM, which decreases from the cluster centre 
to the outskirts by several magnitudes and the inhomoge- 
neous metal distribution. The integration of all quantities 
along the line of sight and the investigated area leads to 
the discrepancy between the true 3D metal mass and the 
metal mass obtained by the mean metallicity of the ICM. 
This is of course an extreme example, but it demonstrates 
how strong the mean metallicity of a whole cluster depends 
on the actual metal distribution. In the case of a constant 
metallicity over the whole ICM, the metal mass obtained 
by the mean metallicity in a circle of 1 Mpc radius around 
the cluster centre is of course exactly the same as the true 
3D metal mass. 

5.2. The mean metallicity as a measure for the true metal 
mass 

In order to test, whether the mean metallicity in a defined 
region around the cluster centre represents a good mea- 
sure for the true metal mass within this area or not, we 
performed several tests with the model clusters. The metal 
content in the model galaxy clusters were calculated from 
an enrichment process investigation, where galactic winds 
and ram-pressure stripping acts simultaneously to enrich 
the ICM. The simulation starts at a redshift of ^^20 and 
results in a metal distribution within the ICM as shown in 
Fig. [71 For more details on the enrichment processes and 
the metal distribution see Kapferer et al. (2007). At red- 
shift z = we extract the mean metallicity within circles 
of different radii around the cluster centre (100, 200, 500 
and 800 kpc) and extract from the 3D simulation the mean 
ICM density within the corresponding spheres. As we know 
exactly the amount and distribution of heavy elements, we 
can compare these two quantities, i.e. the true heavy ele- 
ment mass and the X-ray measured heavy element mass. In 
Fig. [HI the spatial distribution of the density, temperature 
and metallicity of model cluster A is presented. Whereas the 
density does not show significant substructures, the metal- 
licity is inhomogeneous. 

In Table [H the mean metallicity, true metal and X-ray 
measured metal mass for different circles around the clus- 
ter centre for model cluster A and B are given. We find, 
that the X-ray measured metal mass always results in lower 
heavy element masses compared to the true metal mass in 
the same region of interest. If the radius is smaller, the dis- 
crepancy can be up to three times more metals in the ICM, 
than measured by X-ray spectra. By increasing the radius, 
i.e. averaging over much larger volumes, the discrepancy 
becomes smaller. At a radius at 800 kpc the X-ray mea- 
sured heavy element mass gives 20% to 50% lower masses 
than the true metal mass. In Fig. [9l two metal- mass profiles 
for model cluster A are shown. The metal mass increases 
in general with increasing radius, because the volume of 
the radial shell increases with the square of the distance 
to the cluster centre. Although the metallicity decreases 
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Fig. 7. Metallicity profiles for model clusters A and B. 



at larger radii, the larger volume results in a higher metal 
mass. The metal mass X-ray profile (dashed line) is consid- 
erably lower then the true metal mass profile. The discrep- 
ancy arises by averaging of the metallicity in annuli at all 
radii, which wipes out very prominent metallicity peaks all 
over the cluster. In Fig. [8] the distribution of the density, 
temperature and metallicity of the ICM in model cluster A 
at redshift is shown. The high metallicity regions in the 
outskirts are clearly visible in the 3D metal distribution. 
The inhomogeneous metal distribution in galaxy clusters 
leads to a systematic underestimate of the metal mass by 
X-ray observations. For example the ratio of the integrated 
true metal-mass profile to the integrated X-ray metallicity 
obtained metal-mass profile in an 200 kpc radius is 2.70. 
Investigations of the amount of metals residing in galaxies 
in galaxy clusters and in the ICM (e.g. Renzini et al. 1993) 
do always underestimate the true metal mass within the 
ICM. These works concluded, that on average the amount 
of heavy elements within galaxies and the ICM is compa- 
rable, would get a shift towards the conclusion, that even 
more metal mass resides within the ICM as in the galaxies. 

6. Discussion and Conclusions 

We investigate the quality of the metallicities obtained 
by X-ray observation of the ICM as a measure for the 
true metal mass of the ICM. In addition we test if X-ray 
weighted metal maps of galaxy cluster simulations are com- 
parable to synthetic X-ray spectra composed of thermal 
bremsstrahlung and line radiation. The results of our work 
are as follows: 

- If the ICM temperature is above T > 3 x lO'^ K X-ray 
weighted metal maps results nearly in the same metallic- 
ities as obtained from synthetic X-ray spectra. The dif- 
ference is on average less than 7% in areas with r < 750 
kpc around the centre of galaxy clusters. 

— We test how much the distribution of metals within a 
sphere of r= 1 Mpc around the cluster centre changes 
the X-ray measured mean metallicity of this area. We 
found that the X-ray obtained mean metallicity can 
change by a factor of 10 in the extreme case of a dis- 
tinct metal blob residing in primordial ICM, depending 
on the distance from the cluster centre. 
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Fig. 9. Metal mass profile per radial bin for model cluster 
A. Solid line: true metal mass within a three dimensional 
radial bin. Dashed line: metal mass obtained by multiplying 
the metallicity profile taken from the X-ray weighted metal 
map (eqn. 4) with the three dimensional gas mass within 
the same bin (bin size 80 kpc). 



— By using realistic metal distributions in simulated 
galaxy clusters we test the quality of the X-ray obtained 
overall metallicity as a measure for the metal mass. We 
found that depending on the size of the investigated 
area around the cluster centre the true metal mass can 
be three times higher than the metal mass obtained by 
X-ray observations. If the radius increases the discrep- 
ancy gets smaller, but in all cases the true metal mass 
is higher than the metal mass obtained by X-ray obser- 
vations. 

The discrepancies of the true metal mass and the X-ray 
obtained metal mass within the ICM in galaxy clusters can 
be understood in terms of averaging. The observer mea- 
sures integrated quantities along the line of sight, includ- 
ing dense regions in the centre and less dense regions in 
the outskirts. In order to be able to get reasonable X-ray 
spectra of the ICM, observers must always collect X-ray 
photons over an area. As the metallicity is not constant 
throughout the galaxy cluster, i.e. the ratio of primordial 
ICM to heavy elements is not constant, the integration of 
thermal bremsstrahlung and line radiation results in dif- 
ferent ratios, leading to too low metal masses. Applying 
multi-temperature instead of single-temperature models for 
the ICM would alleviate the problem partially as shown by 
Buote & Canizares (1994). 

Taking the results of this work into account, it is obvious 
that estimates from X-ray observations on the heavy ele- 
ments (e.g. iron) content in the ICM lead to masses, which 
are too low. 
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True metal mass [M©] 
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C.3 to C.4 
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1.07 X 10^'' 


8.40 X 10" 
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A 


500 


4.50 X 10^ 


2.63 X lO'' 
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A 


200 


6.05 X 10® 


2.24 X 10* 


0.370 


2.70 


A 


100 
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3.10 X lO'' 
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3.13 


B 
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6.70 X 10'" 


0.94 
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B 
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0.95 
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B 
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B 
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0.97 
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